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ABSTRACT 


We have measured the mean isotopic composition of 
even-charge cosmic ray elements with 14 < Z < 26 near 
0.8 GeV/N using a balloon-borne ionization-chamber/ 
Cerenkov-counter detector system. The experimental 
method makes use of the geomagnetic field as a magnetic 
spectrometer. Results indicate that the most abundant 
isotopes at the cosmic ray source are Si-28, S-32, and 
Ca-40, like the solar system; but Fe-54, unlike the solar 
system. 
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We report here the first results on the isotopic com- 
position of cosmic ray elements of even charge (Z) in the 
interval 14 < Z < 26 for energy near 0.8 GeV/Nucleon. 

Previous workers^"^ have published results on the chemical 
composition in this charge and energy region, however this 
is the first report on the mean atomic weight (A) of these 
elements. Isotopic measurements of cosmic rays, in particular 
that of Fe, have great astrophysi cal significance since they 
have a direct bearing on the nucl eosyntheti c processes which 

create the heavy elements. The technique employed here 

7-8 

was first described by Lund et al . We make use of the 

geomagnetic rigidity cutoff and a detector which measures 

charge and velocity of individual nuclei to determine the 

relative mass to charge (A/Z) ratio for different elements. 

The detector system (Figure 1) is comprised of pulsed 

ionization chambers, a Lucite Cerenkov counter, and a 

scintillation counter hodoscope; the geometry factor is 
2 

0.9 m ster. The system was flown on a balloon launched 

from Southeastern Missouri on 13 September 1972 which 

2 

floated near 3.2 g/cm residual atmosphere. Results on 
chemical composition and energy spectra derived from this 
flight will be published at a later date (see also Ref. 6 
and Ref. 9). The data presented in this paper were gathered 
during a period of nine hours of flight when the balloon 
remained very near a line of constant geomagnetic cutoff 
(~3.0 GV vertical cutoff). This analysis includes about 
4500 events due to Fe . 
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The important character! s ti cs - of the detector system 
for the isotope study are its good charge resolution, good 
velocity resolution, and large geometry factor. The charge 
resolution varies from a = 0.32 charge units at Si to 0.34 
charge units at Fe . Due to the relatively low abundance 
of the odd-Z elements the contamination of the even-Z 
elements is quite small ( 2 % at Fe to 11% at Ar) and has 
been taken into account. 

The Cerenkov counter areal non-uniformities have been 
measured in the laboratory so the velocity resolution is 
limited only by photo-electron statistics. For the range 
of elements and energies with which we are concerned the 
momentum-per-nucl eon resolution is sufficiently narrow 
that it has a negligible effect on our isotope results. 

The isotope analysis procedes as follows. The rigidity, 
p, of a nucleus is related to its momentum-per-nucl eon , p, 
by 

P = (A/Z)p . (1) 


We fit the integral momentum-per-nucleon spectrum for the 
ith element to a power law over a range of p near the cut- 
off. 
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If the total observed flux of the ith element is c , then 
we can define the effective cutoff rigidity, p^ , using 



Eqs . (1 ) and (2) , by 


0. = a. (Z./A. p.) 1 

ic l v i ' i M i ' 
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Combining Eqs. (2) and (3) we can write 
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for any p provided Eq . (2) holds for all allowable momenta 
less than p. Combining Eq. (4) for element i with the 
similar equation for element j we find 
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where e = l/y i - 1 /T j - Although is still not well known 
the ratio p./p. is identically unity if y. = y.. In' fact 

s) * 1 J 

we find the spectra of Fe and Fe secondaries to be somewhat 
flatter than those of Si and S . We find y $ ^ = -1.50 and 
y Fe =-1.43 in the momentum-per-nucleon region near geo- 
magnetic cutoff. Using these y's, a penumbral width of 
15% (typical for cutoffs near 3.0 6V), and a worst case 
penumbral structure we find P-pg/Pg-j differing from unity 
by less than 5 x 10”^. Thus by comparing two elements, 
measured at the same location by the same detector, we 
essentially eliminate the p^ in Eq. (4). Equation (5) 
therefore permits us to determine relative values of A/Z 
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for two elements from the observed total fluxes (J ) above 

c 

cutoff and the fluxes (Op) above some arbitrarily chosen 

momentum-per-nucleon (p). The value of p must be chosen 

high enough to insure Jp is unaffected by the geomagnetic 

cutoff. Since we find the Si spectrum to fit the power law 

9 

down to at least 1.53 GeV/c, we select this as the value 

of p used for determination of He have also analyzed 

the data using larger values of p and find similar results 

but with larger statistical errors, (if the flux of one 

element consists of several isotopes, the A,, in Eq. (5) 

Y k -Y k 

should be replaced by <A^ > where the brackets denote 
the average over the isotopes and the same y^ is assumed 
to hold for all isotopes. However, the precision of our 
results does not warrant this distinction, and we treat 
A^ simply as the mean atomic weight.) 

In Table 1 we present results of the analysis. Column 
1 shows the relative values of A/Z with respect to that of 
Si, errors include statistical uncertainties on the J's 
and the experimental uncertainty on the e. Column 2 shows 
the same data as column 1 corrected for contamination due 
to atmospheric secondaries. This correction is due to the 
small fraction of nuclei (ranging from 6 % of the S to 15% 
of the Cr) observed at the detector which are the results 
of interactions in the atmosphere. These nuclei passed 
through the geomagnetic field as one element, with a 
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specific A/Z, interacted, and were detected as secondaries 
w i th a 1 owe r Z . 

Column 3 shows the resultant values of the mean A 
calculated from Column 2 vnth arbitrary assumption of a 
mean A of 28.5 for Si. Thus far we have determined only 
relative values of A but we note that the relative values 
of A for Si and Fe are not consistent with the solar system 
values.^ If cosmic ray Si is primarily Si-28 then Fe must 
have a large portion of Fe-54. 

In order to determine absolute value of A/Z we need 
prior knowledge of the mean A for some element or elements. 
We use the fact that the observed Ar , Ti, and Cr at earth 
are almost entirely spallation products of Fe on the inter- 
stellar gas.l We have calculated the mean A expected for 
each of these spallation products and shown that the result 
is insensitive to the isotopic composition of the source. 
This result was derived from a propagation calculation 
consisting of solving the set of coupled equations 

N j /X + N j /X j = Q j + p’i /X i )p ij (7) 

where Q. and N. refer to the source strength and observed 

J J 

abundance, respecti vely, of a given isotope of a given element 
p - j are the fragmentation parameters derived from the semi- 
empirical cross-section formulae of Silberberg and Tsao,^ 


assuming unstable fragments with t-j^ < 10® year decay to 

the stable daughter. The are the interaction mean free 
1 3 

paths; X is the mean leakage length for escape from the 
galaxy. By assuming that the elements with 21 < Z < 25 
have a negligible abundance at the source we find that a 
leakage length of 4.5 g/cm allows us to account for the 
observed abundances of each of these secondaries. Using 
this leakage length we can account for all the observed 
elemental abundances for (14 < Z < 26) with a source con- 
sisting only of Si, S, Ca, and Fe in the ratio of 0.97:0.18 
0.15:1.00. (As has previously been noted^’^ there are 

differences between the chemical composition of the cosmic 

1 2 

ray source and the solar system. The solar system values 
are given in Column 6 for comparison.) Results of chemical 
composition from the propagation calculation are shown in 
Column 5. Column 4 shows the observed chemical abundances 
extrapolated to the top of the atmosphere. 

Since Ar, Ti and Cr. are negligible at the source we 
can determine a theoretical value for the mean A of these 
elements at earth from this propagation calculation. We 
find the mean A of these secondaries varies only by 0.3 AMU 
for Cr and 0.1 AMU for Ti and Ar as we vary the source Fe 
from pure Fe-54 to pure Fe-56 and the source Ca from pure 
Ca-40 to pure Ca-48. Our best theoretical estimates for 
these mean A are; Ar-37.3, Ti-47.1, Cr-51.4. Assuming the 
mean A of Ar is 37.3, and using our measured fluxes with 
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Eq. (5), we calculate the mean A of Fe to be 54.3 ± 0.7. 
Similarly, using Ti and Cr as references we find 54.2 ±0.7 and 
55.2 ± 0.7 respectively, in agreement within the statistical 
errors. From these three we derive our best estimate of 
the mean A of Fe at earth to be 54.6 ± 0.5. Similarly using 
the theoretical values of Ar, Ti , and Cr as references, our 
best estimates for the mean A of Si , S and Ca at earth are 
28.5 ± 0.5, 32.1 ± 0.6 and 41.8 ± 0.5. These results are 
summarized in Column 7. 

We have used the same propagation model to investigate 
the source isotopic composition of Si, S, Ca, and Fe. The 
mean values of A at the source to give the observed values 
of A at earth are shown in Column 8. 

The above analysis assumed that all isotopes of a given 
element have the same y. Calcium presents a special problem 
in this regard. We observe y^. = -1.47, midway between that 

of Si (-1.50) and that of Fe and the Fe secondaries (-1.43). 9 
Since the Ca at earth is about 60% primary and 40% secondary 
it seems likely that the observed y^ a can be explained by 
assuming the primary Ca has y = -1.50 and the secondary Ca 
has y = -1.43. Using these assumptions we can account for 
the spectral differences in the analysis. In this particular 
case, since the y's are not very different, the uncertainties 
introduced by this effect are small. 

For the sake of comparison we show in Column 9 mean 

1 2 

atomic masses in the solar system. Note that Si, S and 
Ca at the cosmic ray source are mostly alpha-particle nuclei. 
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similar to the solar system. However, while the solar system 
Fe is overwhelmingly Fe-56,a very significant fraction of the 
cosmic-ray Fe is Fe-54. The implications of this result have 
serious ramifications with respect to the properties of cosmic- 
ray sources and cosmic-ray nucleosynthesis which warrant 
further theoretical and experimental investigation. 
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Alfred P. Sloan Foundation as a Sloan Research Fellow. 

We are indebted to Mr. J. W. Epstein and Mr. A. J. Taylor 
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.and the performance of the balloon flight, and to Mr. A. 

Frank for his assistance in data reduction. The balloon 
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the Office of Naval Research and carried out by the flight 
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F i g u re 1 
Table 1 


A 

CAPTIONS 


Detector System 


Chemical and isotopic composition of cosmic rays 
with 14 < Z < 26 at 800 MeV/N. 
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(A/Z)/(A/Z)Si 
at 3.2 g/cm 2 

(A/ZKA/Z)Si 
Top of Atmos. 

1 

A Relative 
to Si -28. 5 

Observed Chemical 
Composition (Top \ 
of Atmosphere) 

Propagation 
Results X = 
4.5 g/cm 2 

Solar System 

Chemical 

Abundance 

Cosmic Ray 
A at Earth 

Source A 

Solar 
System A 

1 

14 

1.000* 

1.000* 

* 

28.5 

- 1 . 37±0.04 

1.39 

1.12 

28.5±0.5 

28.4±0.5 

28.27 

15 

- 

- 

- 

0.048±0.019 

1 

0.057 

0.014 

- 

- 

- 

16 

0.989±0.01 1 

0.987±0.011 

32.2±0.4 

0.322±0.017 

0.332 

0.569 

32.1±0.6 

31.7±0.9 

32.27 

17 

- 


- 

0.052±0.014 

0.051 

0.002 

- 

- 

- 

18 

1.025±0.015 

1.024±0.015 

37.6±0.6 

0.1 44± 0.015 

j 

0.117 

0.256 

* 

37.3 

- 

36.59 

19 

- 



0.062±0.022 

0.068 

0.004 

: 

- 

- 

i 

| 

20 

1 . 026±0 . 015 

i 1.026±0.015 

41.&0.6 

0.275±0.016 

1 

0.284 

1 

0.082 

41. 8± 0.5 

40.9±1 .0 

| 41.54 

j 

21 

- 

- 

- 

0.059±0.015 

0.039 

0.00004 


- 

j 

22 . 

1 . 056±0.017 

1 . 06Q± 0.017 
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0.003 

* 

47.1 

- 

! 47.20 

! 

23 

- 

- 

- 
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0.001 


- 

I 

[ 

24 

1 .039±0. 016 

1 . 040± 0.016 

50.8+0.8 

0.179±0.013 

0.182 

0.014 

* 

51.4 

- 

j 52.62 
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- 

- 
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| 0.107 '' 

0.010 

- 

- 

j 
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1.032±0.012 

i 
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I 1.000 1 
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1.000 
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